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We study, within the Schwinger-boson approach, the ground-state structure of two Heisenberg 
5— ( , antiferromagnets on the triangular lattice: the J1—J2 model, which includes a next-nearest-neighbor 

coupling J2, and the spatially-anisotropic Ji~J[ model, in which the nearest-neighbor coupling 
takes a different value J[ along one of the bond directions. The motivations for the study of these 
systems range from general theoretical questions concerning frustrated quantum spin models to the 
concrete description of the insulating phase of some layered molecular crystals. For both models, 
the inclusion of one-loop corrections to saddle-point results leads to the prediction of nonmagnetic 
phases for particular values of the parameters J1/J2 and Ji/Ji . In the case of the Ji — J2 model we 
' shed light on the existence of such disordered quantum state, a question which is controversial in the 

literature. For the Ji — J[ model our results for the ground-state energy, quantum renormalization 
of the pitch in the spiral phase, and the location of the nonmagnetic phases, nicely agree with series 
expansions predictions. 
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The triangular-lattice Heisenberg antiferromagnet (TLHA) has played a fundamental role in the understanding of 
frustrated quantum spin systems. In particular, the existence of a nonmagnetic ground state in the S = 1/2 system 
[ has been stronglypdebated in the literature, although there is now a widespread conviction that it displays the classical 
120° spiral order.Elia In the last years, this particular problem was linked to more general questions concerning spin- 
O liquid states and their possible connections to high- Tc superconductivity. In this context, the most studied model has 
I been tbe square-lattice Heisenberg antiferromagnet with first- and second-neighbor interactions, the so-called Ji — J2 

model.H The introduction of frustrating second-neighbor interactions leads, at intermediate values of the couplings, 
y—i to the existence of a disordered spin- liquid state which intervenes between the Neel and "striped " (ferromagnetic in 
^ ' one direction, antiferromagnetic in the other) magnetic orders. This situation is by now fairly well established by 
a variety of numerical and analytical methodsJju On the other hand- an extension of the TLHA including second- 
neighbor interactions, the J1 — J2 TLHA, has also been considered.Q^llil This extension is a natural development after 
the thorough investigation of the model on the square lattice, but the existence or not of an intepHiediate spin-liquid 
i phase in this case is not so clear, with most works favoring the non-existence of such a state.Q^tJ This is somehow 
Q\ ' paradoxical, since the TLHA was considered for many years the best candidate to have a spin-liquid ground state 
0^ ■ just on the basis of the frustrating lattice topology. One would expect that the introduction of additional frustration 
' through the J2 interaction should contribute to melt the already weak order of the nearest-neighbor model. 
, The study of the J1 — J2 TLHA was mostly driven laLnuxely theoretical motivations. On the other hand, recent 
^ • experimental results have produced a surge of interestll3^llj in the ground-state properties of the nearest-neighbor 
TLHApwith spatially anisotropic couplings: Ji along two bond directions and J[ in the third one. It has been 
arguedEj that this model should describe the magnetic phases of the quasi-two-dimensional organic superconductors 
Q ■ AC— (BEDT-TTF)2X. For these layered molecular crystals the relevant values of J'l/Ji are around 0.3 — 1.0, with the 
^ , ratio varying with the anion X and with uniaxial stress along the layer diagonal. 

ILjl ■ We have previously performed a study of both the -^i and Ji—J'i TLHA using the rotational invariant Schwinger- 
boson approach (SBA) in a mean-field approximation.EI At that time, our motivations for studying the latter model 
was the natural interpolation that it provides between the nearest- neighbor TLHA {J[ = Ji) and the square-lattice 
] antiferromagnet (J( = 0). We obtained a good agreement between the ground-state energy predicted in our approach 
. - • and exact numerical values on finite lattices. Furthermore, no indicationj-pf a disordered state was found for the values 
of J2/J1 and J[/Ji of interesLi-'in this work we extend the calculations irJj to include one- loop corrections to the mean- 
field picture. We have shownaO that these corrections bring the saddle-point results in line with exact diagonalization 
values on finite clusters, which lends support to the SBA predictions for the thermodynamic limit. In particular, for 
the square-lattice J1 — J2 model we found that there is a quantum nonmagnetic phase for 0.53 ^ J2/J1 ^ 0.64. This 
result was obtained by considering the spin stiffness on large lattices and extrapolating to the thermodynamic limit, 
a procedure that avoids the infinite-lattice infrared divergencies associated to Bose condensation. We show here that 
similar behaviors are predicted for both the Ji — J2 and Ji — J[ TLHA. 
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II. THE CALCULATIONAL METHOD 



We consider a general Hamiltonian of the form 



i5]j(r-r')^r-^r', (1) 



^ 2 

r,r' 

where r,r' indicate sites on a triangular lattice. As usual, we write spin operators in terms of Schwinger bosons:lii 
Si = ia^.a.ai, where a! = (aj^a-j^) is a bosonic spinor, a is the vector of Pauli matrices, and there is a boson- 
number restriction 0L\aO-ia = 25* on each site. With this representation the spin-spin interaction can be written 
as Si.Sj —: BjjBij : —AljAij, where we defined the S'C/(2)— invariant order parameters Aij = cra,i„ajg and 

^\j ~ 1 X](T '^L'^jo- — —(J, a = ±1), and the notation .O: indicates the normal order of operator O. Thus, we can 
formulate the partition function for the Hamiltonian (Q) as a functional integral over boson coherent states, which 
allows its evaluation by a saddle-point expansion. Since the theory presents a local C/(lX-symmetry, we use collective 
coordinate methods — as developed in the context of relativistic lattice gauge theoriealZl — to handle the infinitely- 
many zero modes associated to the symmetry breaking in the saddle point. These modes without restoring forces, 
which correspond to local gauge transformations, are eliminated by the exact integration of the collective coordinates. 
Such a program can be carried_out by enforcing in the functional integral measure the so-called background gauge 
condition, or "natural" gauge,EZl by means of the Fadeev-Popoff trick. In this way we restrict the integrations in the 
partition function to field fiuctuations that are orthogonal to the collective coordinates. At T = 0, after carrying out 
the integrations on these genuine fluctuations, we obtain the one-loop correction to the ground-state energy per site, 

E.^-^rd.j:iJ |. (2) 

Here the Fadeev-Popoff determinant AFp(k,u;) — |(^q (k,a;).(/9^(k, w)] , where i^Q(k,Ci;) is the left zero mode of the 

dynamical matrix A^'^'^ in the k — subspace, and A^^^ is the projection of this matrix in the subspace orthogonal 
to the collective coordinates. The dynamical jpatrix ^'^^^ is the Hessian of the effective action as a function of the 
decoupling (Hubbard-Stratonovich) fields (seeQ for details). 

In the ordered phases_af the model the Bose condensate breaks the global SU{2) symmetry and its density gives 
the local magnetization£3 The associated physical Goldstone modes at k = 0, ±Q (Q is the magnetic wavevector) 
lead to serious infrared divergencies of intermediate quantities, which have to be cured by standard renormalization 
prescriptions. In order to avoid these problems we have computed physical quantities (which are free of divergencies) on 
large by finite lattices, and finally extrapolated these values to the thermodynamic limit. We considered clusters with 
the spatial symmetries of the infinite triangular lattice, corresponding to translation vectors Ti — (n + m)ei + me2, 

T2 — nei + {n + m)e2^^ Here ei ^ (a, 0) and 62 = (— ^a, -^a) arc the basic triangular lattice vectors. To fit to 
the cluster shapes the expected (spiral) magnetic orders, and also to allow the calculation of the stiffness, we impose 
arbitrary boundary conditions 5'r-i-Ti = 'Tli^{^i)Sr {i = 1,2) , where TZ-^{^i) is the matrix that rotates an angle $i 
around some axis n (notice that we are using boldface (arrows) for vectors in real (spin) space). It is convenient to 
perform local rotations S^ TZ-^{9r)Sr of angle 9^ = AQ • r, so that with the choice AQ • — the boundary 
conditions become the standard periodic ones S^+t, = S^- Then, we define the (T = 0) stiffness tensor by 



9^i?Gs(Q+AQ) 
pA = 



(3) 

AQ=0 



where Eqs is the ground-state energy per spin and 9i = AQ • (i = 1, 2) are the twist angles along the basis vectors 
e^. In the next two sections we will apply this formalism to the models under consideration. 

III. THE J1-J2 TLHA 

The J1 — J2 TLHA is given by (|l|) with J(r — r') = Ji(J2) for r and r'nearest (next-nearest) neighbor sites, and 
otherwise. For classical spin vectors, when a = J2/J1 < 1/8 the lowest-energy configuration in this model is the 
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commensurate spiral with magnetic wavevector Q = (|^,0); for 1/8 < a < 1 there is a degeneracy between the two- 
sublattice and fom'-sublattice Neel orders. Quantum fluctuations Uft this degeneracy and select, through an "order 
from disorder " phenomenon, the two-sublattice coUinear ground-state with magnetic wavevector Q = (t^-^^^^ 
This scenario was first proposed using spin-wajve theoryli and later confirmed by a study of the thermodynamic-limit 
collapse to the ground state of low lying levelsO The correction (^) to the saddle-point value Eq leads to the ground- 
state energy Eqs = Eq + Ei shown in Fig. 1. This figure contains the result for the infinite lattice and also for a 
finite cluster of 12 sites, which allows a comparison with exact results obtained by numerical diagonalizationJj We see 
that the addition of the Gaussian correction Ei improves the saddle-point value Eq, particularly for the 120° spiral 
phase. Moreover, the departure of the fluctuation-corrected results from the exact values in the range 0.1 ^a^S 0.2 
hints to the possible existence of a disordered phase in this region. In the thermodynamic limit, at saddle-point order 
the theory predicts a first order transition between the two magnetic ground states at some intermediate frustration 
a ~ 0.16, with no intervening disordered phase.EI After the inclusion of the Gaussian fluctuations there is a window 
0.12 ^ a < 0.19 where the stiffness vanishes (see below) and the magnetic order is melted by the combined action 
of quantum fluctuations and frustration. This result should be compared with the linear spin wave results oiu, which 
predict a quantum nonmagnetic phase in the range (0.1,0.14). Notice, however, that within spin-wave theory this 
window closes when corrections to the leading-order calculations are included.^ We stress that the same happens in 
this approach for the Ji — J2 model on the square lattice, where however other methods confirm the existence of a 
nonclassical phase. 

As stated above, the existence or not magnetic long-range order in the thermodynamic limit was investigated by 
considering the spin-stiffness tensor (^).c3 For spins lying on the a;y- plane, it is necessary to consider both the parallel 
stiffness = under a twist around z and the perpendicular stiffness p± = p^^ for twists around an arbitrary versor 
n on this plane. However, on clusters with periodic boundary conditions our approach is rotational invariant and we 
only have access to p = ^{P\\ + 2/0_l). On the other hand, we found that the large-S* Schwinger-boson prediction for 
this quantity is exactly 4/3 smaller than the corresponding classical result. As discussed inlla, to solve these problems 
we considered clusters that fit the magnetic orders in the xy-plane with appropriate twisted boundary conditionsJlju 
In this case, since the rotational invariance is explicitly broken by the boundary conditions, one is able to compute the 
parallel stiffness and, -moreover, the large-5' predictions have the correct behavior (no ad hoc factors are required 
to correct this quantity) .113 Finally, using the values for p obtained from clusters with periodic boundary conditions it 
is possible to determine pj^ = ^ p — ip||. The corresponding results are presented in Fig. 2, both at saddle-point and 
one-loop orders for comparison. Notice that in the 120° spiral phase py softens first than p±, while the behavior is the 
opposite for the colincar state. Since the indirect calculation of p± requires the use of ad hoc factors to renormalize 
p, the determination of the upper limit for the disordered phase might be unreliable. However, based on our previous 
experience with the SBA, we believe the existence of an intermediate nonmagnetic phase can be trusted. 



IV. THE Ji-Ji TLHA 

For the Ji - J[ TLHA, in (|l|) we take J(r - r') = Ji for r - r' = ei (i = 1, 2) and J(r - r') = J[ for r - r' = 63 = 
ei +62. As mentioned above, this model is interesting in viessz. of its connections to the spin degrees of freedom in 
the insulating phase of some layered organic superconductors.!^ Very recent works havoipvestigated its ground-state 
phase diagram and other properties, using spin-wave theorylij and series expansions.Eil For classical vectors, with 
T] = J'l/Ji < 1/2 the lowest-energy configuration is the two-sublattice Neel order discussed in the previous section, 
with the (frustrated) ferromagnetic arrangement of the spins along the weakly-coupled 83 direction. For ry > 1/2 
the preferred spin configuration becomes an incommensurate spiral, with the angle Qciass between neighboring spins 
along the ei,e2 directions given by Qciass = arccos(— ^). 

The inclusion of the quantum nature of spins by means of the SBA produces the results for the ground-state energy 
shown in Fig. 3. In this figure we give the ptnean- field values and the fluctuation-corrected results, and compare them 
with the series expansion predictions fromliJ. We see that the fluctuations produce regions in which the saddle-point 
solutions become unstable. On the other hand, for those values of rj where the magnetic phases considered are stable 
there is a very good agreement between both methods. It is also of interest to compare the quantum renormalization 
of the spiral pitch; in Fig. 4 we plot the classical result for Q given above, and the corresponding angle that minimize 
the quantum ground-state energy. Again, in the region where the magnetic phases are stable the results nicely agree 
with those coming from series expansions. Furthermore, we have checked at mean-field order that the angle that 
minimizes the total energy corresponds also to the minimum of the quasiparticle dispersion relation, a fact that was 
used irO to determine Q. 

To establish the regions without magnetic order we studied again the spin-stiffness behavior. Like for classical 
vectors, the stiffness tensor (^) is diagonal in the (perpendicular) directions 83 and 82—81; the corresponding stiffness 
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along these directions are plotted in Fig. 5. We found a qualitatively different behavior between the saddle-point 
and one-loop results; most notably, in the Neel phase the Gaussian fluctuations soften first the classically stronger 
stiffness in the 63 direction. At mean-field order our calculations indicate a continuous transition between the two 
magnetic phases at 77 ~ 0.62f , and the absence of a magnetic saddle-point solution beyond 77 ~ 2.20. The corrected 
results show the melting of the Neel phase at r; ~ 0.61, still above the classical point ?7ciass — 1/2, and predict a 
disordered quantum phase in the range 0.61 ^ rj 0.96. Furthermore, the incommensurate phase becomes stable 
only in a reduced parameter region 0.96 ^ 77 ^ 1.10. In this case the instability appears as a negative eigenvalue of 

the jpcpjected dynamical matrix A^^^ in (Q). These results are in fair agreement with the series-expansion predictions 
fromEJ, which indicate that the Neel order disappears at ~ 0.65 — 0.7, the disorder region extends from this value 
up to rj ^ 0.9, and there is an incommensurate phase for rj ^ 0.9 with no clear ending point. 



V. CONCLUSIONS 



In conclusion, we have considered, within the SBA, the Gaussian-fluctuation corrections to the spin stiffness of 
the Ji — J2 and Ji — J[ TLHA. For the Ji — J2 model we found that the order-parameter fluctuations weaken the 
stiffness, which is reduced by increasing the frustration until it vanishes leaving a small window 0.12 ^ a ^ 0.19 
where the system has no long-range magnetic order. Like in previous studies, we found that the consideration of 
clusters which require twisted boundary conditions to fit the magnetic orders avoids the use of ad hoc factors to 
correct the Schwinger-boson predictions. This fact points to a subtle interplay between rotational invariance and the 
relaxation of local constraints in this approach. 

In the case of the Ji — J{ TLHA our results indicate a rich phase diagram, with at least two magnetic (Neel and 
incommensurate spiral) phases and two disordered quantum states in the parameter region of interest. One of the 
most notable aspects of our calculations is the quantum renormalization of the magnetic wavevector in the spiral 
phase, which agrees remarkably with the series expansion prediction. We also found that the Neel order extends 
beyond its classical stability point up to a value rj ~ 0.61, where it seems to melt continuously into a purely quantum 
phase. On the contrary, the spiral order is favored in a reduced parameter range, with our results indicating first 
order transitions from this phase to the nonmagnetic states. These first order transitions appear in our calculations 
as local instabilities against the order-parameter fiuctuations. 

Finally, there still remains the difficult task of characterizing the disordered quantum phases of the modeJs under 
consideration. Some attempts in this direction have already been done, for both the Ji — J2B and Ji — J{Eil TLHA, 
but there is not clear understanding of these phases yet. They are usually considered to be of the resonantp*£ilence- 
bond type, and are in general described starting from a regular strong-bond (dimer) covering of the lattice.til These 
investigations can be performed within the formalism developed here, since the SBA does not rely on having magnetic 
order in the system like, for instance, spin-wave theory. However, these studies would require the extension of the 
present calculations to larger magnetic unit cells and the computation of physical quantities able to characterize the 
new phases, a question that is far from trivial. 
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FIG. 1: Ground-state energy per site Eqq of the Ji — J2 TLHA as a function of the frustration parameter a — J2/ Ji- 
Dashed and full lines give the mean-field and fluctuation-corrected results, respectively. The lower curves correspond 
to the energy of a 12-site cluster (the values are shifted by -0.05 for clarity); the upper curves give the thermodynamic- 
liniit results. Dots are exact numerical values fromtl. 

FIG. 2: a) Parallel stiffness py , and b) perpendicular stiffnes p± for the Ji — J2 TLHA as a function of the frustration 
parameter a. Dashed and full lines give the saddle-point and fluctuation-corrected results, respectively. The thin 
vertical lines separate the regions with magnetic order from the middle parameter range where there is a quantum 
disordered phase. 

FIG. 3: Ground-state energy per site Eqs of the Ji — TLHA as a function of — J^j J\. Dashed and full lines 
give rtije mean-field (MF) and fluctuation-corrected (FL) results, respectively. Dots are series expansion predictions 
fromM. 

FIG. 4: Relative angle Q between nearest-neighbor spins as a function of ?/. The thick full line corresponds to 
the lowest-energy configuration of classical spins, the dashed line is the mean-field prediction in the quantum case, 
and the thin full line is the series expansion result fromtil. Open dots are the fluctuation-corrected results in the 
incommensurate phase (shown in more detail in the inset). 

FIG. 5: Spin stiffness py for the Ji—J^ model as a function of 77. Dashed and full lines correspond to saddle-point and 
fluctuation-corrected results, respectively. The thin vertical lines separate the regions with antiferromagnetic (AF) 
and incommensurate (INC) spiral orders from the parameter ranges corresponding to disordered quantum phases. 
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